 

------------------------------------------------------------------------------------------------------------------------------------------------
       Step 1: Compute and print the minimum required probe dimensions
------------------------------------------------------------------------------------------------------------------------------------------------

Since the cross section Th232(p,n3)Pa230 becomes 0 for energies below the 15 Mev threshold,
   Usefull Probe Thickness: 0.165 cm
   Usefull Probe Diameter: 181 um + beam diameter


------------------------------------------------------------------------------------------------------------------------------------------------
       Step 2: Calculate the average cross section of Th232(p,3n)Pa230 over the energy range
------------------------------------------------------------------------------------------------------------------------------------------------

      Diagram:
Cross Section[mbarn],  and Stopping Power[keV/(g.mm2)]
Energy [MeV], XSection [mbarn], -dE/dx [keV/(g.mm2)]
15, 0, 130.158
16.4, 90.5, 122.554
17.6, 149.5, 116.851
18.2, 239.1, 114.149
19.9, 352.5, 107.346
21, 342.9, 103.644
21.5, 318.7, 101.943
22.5, 315.9, 98.591
23.5, 280.7, 95.800
24.3, 217.9, 93.568
25.2, 169.5, 91.147
26.6, 107.2, 87.835
28, 81.8, 84.694
28.6, 80.2, 83.473
31.6, 64.8, 77.801
33.5, 53.1625, 74.649
34, 50.1, 73.879

Average Cross Section of Th232(p,3n)Pa230:
   Usefull Mass Thickness: 1.918 e+00 g/cm2
   <XSection> = 162.63 mbarn


------------------------------------------------------------------------------------------------------------------------------------------------
      Step 3: Calculate the Production Rate of Pa230
------------------------------------------------------------------------------------------------------------------------------------------------

The reaction rate with a beam of 200 uA is:
       ReactionRate  =  I / e * XSectionAv * 1e-24 * MassThickness / AW * Const_Avogadro [Reactions / sec]
   ReactionRate = 1.010 e+12 Reactions/sec
   Proton flux: 1.248 e+15 protons/sec


------------------------------------------------------------------------------------------------------------------------------------------------
      Step 4: Calculate the Production of Pa230 and U230 over 50 hours
------------------------------------------------------------------------------------------------------------------------------------------------

After 50 hours beaming at 200 uA, we get:
   1.745 e+17 atoms Pa230 = 8.047 e+01 GBq     * This is quite good comparing to the 84.1 GBq from the original paper
   6.027 e+14 atoms U230 = 2.324 e+02 MBq

   Diagram:
Pa230 production trough Th232(p,3n)Pa230
Time (hour), Pa230 (Bq), U230 (Bq)
0, 0, 0
5, 8351254085.91113, 2432781.25485408
10, 16633486452.102, 9681858.59404102
15, 24847267551.554, 21673965.3099589
20, 32993163122.5503, 38336671.7321434
25, 41071734227.6415, 59598376.7072991
30, 49083537292.29, 85388299.1616469
35, 57029124143.1945, 115636469.743295
40, 64909042046.2985, 150273722.544781
45, 72723833744.4838, 189231686.905161
50, 80474037494.953, 232442779.290309


------------------------------------------------------------------------------------------------------------------------------------------------
      Step 5: compute the optimal cooling time for a maximum production of U230
------------------------------------------------------------------------------------------------------------------------------------------------

Optimal cooling time for maximal U230 production: 26.88 days
Neglecting the initial U230 amount: 27.41 days


------------------------------------------------------------------------------------------------------------------------------------------------
      Step 6: U230 production by decaying a nuclide mixture during 27 days
------------------------------------------------------------------------------------------------------------------------------------------------

After 27 days cooling time:
<Image>
ba909755-3c15-4ac2-aeb4-d997d1d1f773.gif

91 Pa230
   Activity: 2.844 e+10 Bq
   Mass: 2.356 e-05 g
   Gamma emission rate: 0.000 e+00 keV/s
   EDC Ingestion: 9.200 e-10 Sv/Bq
   Committed effective Dose: 2.617 e+01 Sv
   ALI for Ingestion: 2.174 e+07 Bq
92 U230
   Activity: 2.350 e+09 Bq
   Mass: 2.327 e-06 g
   Gamma emission rate: 2.464 e+09 keV/s
   EDC Ingestion: 5.600 e-08 Sv/Bq
   Committed effective Dose: 1.316 e+02 Sv
   ALI for Ingestion: 3.571 e+05 Bq


------------------------------------------------------------------------------------------------------------------------------------------------
      Step 7: U230 Packaging and Transport
------------------------------------------------------------------------------------------------------------------------------------------------

U230 2.35 GBq, form: 0, sealed: 0

   Number of Sources: 1

Totals for Source 0:
   Mass : 2.327 e-06 g
   Activity: 2.350 e+09 Bq
   Heat: 2.253 e-03 W
   Gamma dose rate at 1 m: 3.235 e-01 uSv/h
   Fissile Uranium ratio: 0.000
   Fissile Plutonium ratio: 0.000
   Container Type: 0
   Package: Expected package

------------------------------------------------------------------------------------------------------------------------------------------------
